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A new Green's function representation is employed in a microscopic derivation of a Ginzburg- 
Landau theory of strongly type superconductivity at high magnetic fields. An exact analytical, 
physically transparent expression for the quartic term in the corresponding order parameter expan- 
sion is presented. The resulting expression reveals singular non-local contributions to the supercon- 
ducting (SC) free energy, associated with highly coherent cyclotron motions of the paired electrons 
near the Fermi surface, which are strongly coupled to the vortex lattice. A major part of these con- 
tributions arises from incoherent scattering by the spatially averaged pair-potential, which is purely 
harmonic in the de Haas van Alphen frequency. However, coherent scatterings by the ordered vortex 
lattice generate, at low temperatures, large erratically oscillating (i.e. paramagnetic-diamagnetic) 
contribution to the SC free energy as a function of the magnetic field. Vortex lattice disorder, which 
tends to suppress this oscillatory component, is found to preserve the singular harmonic part of the 
SC free energy. 

PACS numbers: 74.25.Ha, 74.25.Uv 



I. INTRODUCTION 

Revealing the mechanism in which a clean, ordered su- 
perconducting (SC) material at very low temperatures re- 
sponds to the application of an external magnetic field is 
of crucial importance for understanding the phenomenon 
of superconductivity at its most fundamental level. Sur- 
prisingly, as compared to high-field superconductivity 
under "non-ideal" conditions, i.e. in dirty or disordered 
materials at relatively high temperatures, the "ideal" SC 
state at low temperatures under high magnetic fields is 
currently not well understood even within the conven- 
tional BCS theoryi^. Experimental detection of such 
high magnetic field states is currently restricted mainly 
to magneto-quantum oscillations techniques^, which pro- 
vide researchers with detailed quantum information of 
the system investigated, but are difficult to analyze by 
the standard Fourier transform techniques due to the 
highly restricted range of magnetic fields available be- 
low the SC transition. There have been many reports on 
observation of de Haas van Alphen (dHvA) or Shubnikov- 
de Haas (SdH) oscillations in the SC states of strongly 
type-H superconductors, showing occasionally additional 
damping of the signal in the SC state with respect to the 
normal state signal (a partial list of references includes 
Refs4"— ). However, their interpretation has not lead so 
far to any kind of consensus regarding the infiuence of the 
vortex matter on the oscillatory signal, partly because of 
the difficulties in the data analysis, and partly due to 
lack of a consistent theory with a quantitatively predic- 
tive power—. The mean-field theories based on a detailed 
exposition of the quasi-particle excitations obtained by 
solving the corresponding Bogoliubov-de Gennes (BdG) 
equations for an ordered vortex lattice^i^— , provide in- 
sight into fine features of the Landau band structure, but 
lose their transparency very quickly and become heavily 



numerical at early stages of their application to any ob- 
servable quantity. On the other hand, a simple formula 
for the additional dampin g^^'^^ , used frequently in the lit- 
erature for a manageable interpretation of experimental 
data, has been shown to be limited to situations of ran- 
dom vortex distributions^, and to the influence of the SC 
order parameter on the quasi particle relaxation rate^^, 
neglecting important contributions to the oscillatory SC 
free energy. 

Attempting to compare the results of the different the- 
oretical approaches leads to great confusion. In the field 
range near where the SC order parameter is small 
and the leading (quadratic) term in the order parame- 
ter expansion of the SC free energy should be a good 
approximation^, the results of all numerical simulation of 
the BdG cquations^rJ^ deviate markedly from this uni- 
versal limiting value (see e.g. Fig. 8 in RefJi^).Thc sit- 
uation in the low magnetic fields region well bellow Hc2 
is even worse. Norman and MacDonald (NM) in their 
numerical simulations of the BdG equations reportedi^ 
that the harmonic of the Fourier transform of the cal- 
culated magnetization in a finite interval of this regime 
varies erratically in sign and magnitude (see Fig. 5 there) 
with no pattern that they could discern. Similar cal- 
culations carried out by Yasui and Kitapi^ have resulted 
in an additional damping rate which deviates drastically 
from both NM result and Maki's formula, showing er- 
ratically oscillating patterns of the corresponding Dingle 
plot, which seems to be indicative of some fundamental 
problem of numerical instability. 

It is therefore desirable to derive an exact analytical ex- 
pression for the SC thermodynamic potential which will 
enables one to identify the origin of this erratic behavior 
and carry out the calculation in a controlled fashion. In 
this paper we present such an exact analytical expression 
within the framework of the Gorkov-Ginzburg-Landau 
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expansion of the free energy in the SC order parame- 
ter (or pair potential) up to fourth order. It is indeed 
found that highly correlated portions of the electronic 
cyclotron orbits involved in pairing have dominant con- 
tributions to the quartic and higher order terms of the SC 
free energy. The corresponding, singularly divergent dis- 
tribution of the free energy obtained within an effective 
temporal representation, results in equally shared contri- 
butions from the spatially uniform component of the SC 
pair-potential and from its periodically modulated com- 
ponent associated with the vortex lattice. 

The resulting expression of the free energy consists of 
two types: terms harmonic in the dHvA frequency and 
associated with the Landau level structure of the quasi 
particles, and erratically oscillating terms as functions of 
the magnetic field, associated with coherent scattering 
channels of the quasi particles by the vortex lattice. The 
latter are strongly enhanced when the length of a recip- 
rocal vortex lattice vector coincides, or nearly coincides 
with the Fermi sphere diameter. 

The resulting analytical expression also enables one 
to study the effect of disorder in the vortex lattice on 
the oscillatory free energy. In the white noise limit of 
the self-consistent Born approximation (SCBA)^^, the er- 
ratic oscillations associated with the coherent scattering 
by the vortex lattice arc suppressed, yielding a result 
consistent with Maki's formula^. However, deviations 
from the SCBA should be carefully examined in light of 
the recent high-field low-temperature ^SR measurements 
in the vortex-glass (peak-effect) region of borocarbide 
superconductors^^, which have shown strong correlation 
betwen the enhanced additional damping of dHvA oscil- 
lations observed in the peak-effect region with enhanced 
vortex lattice disorder in this region^. 

II. GENERAL FORMULATION 

We consider a 2D strongly type-II (neglecting the effect 
of SC screening currents) superconductor in a perpendic- 
ular uniform magnetic field H = (0, 0, H). Generalization 
to isotropic 3D systems is rather straightforward. It is 
assumed that the superconductor can be described by 
means of BCS-hamiltonian density for the usual singlet 
s-wave electron pairing 



a general vortex lattice state, 

written in terms of a discrete set of ground-state Landau 
orbitals: 

Mx,y) = e"^^e-'^"'+*«"^-(^+«"/2)= (2) 

n 

n 

where z ~ x + iy and qn = —n = qon with the lattice 
spacing Qx along the x - axis and the angular parameter 
9 which determines the point symmetry of the vortex lat- 
tice. It is easy to see that for a general (rhombic) vortex- 
lattice geometry, determined by the angular parameter 

9 , = IT / ^/f — (pj-K^ . For the Abrikosov triangular 

lattice: 6 =- t:/2 and go = = 3^/-*V2^. 

We use the symmetric gauge with vector potential 
A=i[Hxr] = i7J(— ?/,x,0) and dimensionless space 
coordinates measured in units of the electronic magnetic 
length, aH = ^/ch/eH . The amplitude of the order pa- 
rameter, Aq = 5"^ / (pTi |A(ri)|^ , where S = ttN is the 
area of the 2D superconductor and A*" is the number of 
vortices, is treated as a variational parameter for mini- 
mizing D,sc (A). 

The thermodynamic potential, flgc (^)j can be written 
as a Taylor expansion in the SC order parameter—: 

n^ci^o) = s- — + V ^^i^2„ (Ao) , 

(huc) gscs n=l " 

^2n = n^^^ j £ {V}f2n{{v})k2n{{v}) , 

where: 

k2n{{v}) = (2^;iV™)'"E^S(ri,r2,o..) X (3) 
G'o(r2, ra, cj,,)...G5(r2„_i, r2„, uj^)Go{v2n,ri-,^v) 



Hbcs = A*(r)V(r)V(r)+A(r)V*(r)V*(r) 
^ |A(r)|^ 



9BCS 



(1) 



where ip (r) is the single electron field operator and gscs 
is BCS coupling constant (electron spin is neglected for 
the sake of simplicity). Within mean- field approximation 
the order parameter, A (r) = gscs (V' (r) tp (r)), should 
be determined self-consistently by minimizing the super- 
conducting thermodynamic potential, r2sc(A). We do 
not consider the problem in a fully self-consistent man- 
ner, assuming that the order parameter is described by 



r2«({r}) = .g*(ri,r2)5(r2,r3)...g*(r2„_i,r2„) x (4) 
5 (r2„ , r 1 ) (r 1 ) <Po (i"2 ) • ■ •¥'0 (r 2„- 1 ) (i"2n ) 

Here we use the normal state single electron Green's 
function in the uniform magnetic field, which is 
given by the well known expression, G'(r2, ri, cjjy) = 
.g(r2,ri)Go(r2,ri,w,,), where: 

Go(r2,ri,cj^) = Gq {p,uj^) = ^^^^ ^ ■ 



Pf — £n + ihjJu 



(5) 
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p = r2 — ri , and g(r2,ri) = e-i["'2>^"'i] is the usual 
gauge factor. Also note that in the above equations 
Lo^ — nksT {2^ + 1) /h, 1/ = 0,±1,... is the Matsubara 
frequency at temperature T, is the chemical poten- 
tial, and Sn = fi^jJc {n + 1/2) , n = 0, 1, 2, ...is a Landau 
level energy with lOc = eH/mc the cyclotron frequency. 

It will be very helpful to use an equivalent represen- 
tation of the translational invariant part of the Green's 
fmiction for cj^ > Oi.e.: 



Go (P) 



1 



27rfixj, 



c Jo 



dre 



dT 



2'KfUjJc Jo 



+e- 



1 -e- 



1 



(6) 



which can be easily derived from Eql5]by using the inte- 
gral representation of {^ip — + itiWy)~^ and the gen- 
erating function of the Laguerre polynomials. The re- 
sulting expression is written in terms of the following 
dimensionless quantities: a = 1 — 
(/ii? — £„) = hucinp — n) and io„ 



1 l+e 



4 l-e-»^ ' 

hw^Wi, . Note 

also that for vOi, < 0, r in EqlH] should be replaced with 
— r , yielding the complex conjugate of the expression for 
Wv > 0. 

Exploiting the integral representation, Eq|6l we can 
rewrite the electronic kernel K2n{{j^}) in the form: 



/>oo 

^2„({r})-n/ dT,e-^''^^ 
„ JO 



1 



a 



(-1) 



a, = 1 - e*' 



1 1 
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This representation of D,sc (Ao) has an obvious advan- 
tage over the original expression: all space integrals are 
of Gaussian forms and, therefore, can be calculated an- 
alytically. Unfortunately, gauge factors mix all electron 
coordinates so that the calculation of the higher order 
terms is not trivial. 



III. QUARTIC VERSUS QUADRATIC TERMS: 
EFFECT OF THE VORTEX LATTICE 



A. The quadratic term 

The second order term have been calculated long ago. 
We repeat the calculation to illustrate the advantage of 
using the Green's function in the special representation, 
EqlH 

The vertex part in Eq|4] can be written as 

r2(ri,r2) = .g*(ri,r2)5(r2,ri)v3o(ri)</'S(i"2) 



(2) 

where Qirh = i {xiy2 - yiX2) + ixiyi - ix2y2 + 
iqnXi - {yi + g„/2)^ - iqm.X2 - (2/2 + gm/2)^. Not- 
ing that the dependence on the center of mass 
coordinates, R = '^^'^^^ , appears only in the ver- 
tex part, one can extract this dependence from 
(2) 

Cnm, ending with two integrals over the center of 
mass: J dRxCx:p[i{qnRx - QmRx)] a^NxSnm and 
/ dR'y exp {-2R'^) = Here is the number of 

vortices along the x-direction and R!y is a shifted Ry 

coordinate. The remaining function, C^m ^\p^ 
does not depend on n, so that summation over n gives 
the number of Landau orbitals along the y-axis, TVy. 
Consequently, the quadratic term can be written as 



inF(T2-Ti)-ro„(Ti+T2) 



00 /'OO 



dT\dT2 



(7) 
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aiQ2 



d^pexp ( -/-tip^ - \i2P^ - T,P^ 



or, after integrating over relative coordinates, p, as 

dTxdT2 (8) 



^2y = a^J'^Nnf^ 



00 /'OO 











m_F(T2-Ti)-ro„(Ti+r2) _ 



1 



ai + a2 



where = NxNy is a number of vortices in the system. 

The dominant contributions to the r-integrals origi- 
nates in the poles of the integrand where ai^2 0, 



' — exp [-fijpj) , namely at t. 



2njTr, rij 



0, ±1, .., where the first expo- 
nent inp {t2 — Ti) is equal to 2i'iTnpn, n = 0, ±1, .., cor- 
responding to exact harmonics of the dHvA frequency 
F = npH. We therefore conclude that the quadratic 
term is dominated by harmonics of the dHvA frequency 
which implies that to leading order in the GGL expansion 
the Landau levels structure is not distorted by the vor- 
tex lattice. This result is consistent with the well known 
property of the quadratic term to be independent of the 
vortex lattice structure. 

Considering the first harmonic for the sake of illus- 
tration, we shift T2 — >■ 27r -t- T2 and expand ai + 0:2 
in Ti and T2 for ri ^ 1 and \t2\ ^1: ai + 0,2 — 
i {t2 - Ti) + \ (ri + T2) . Here the term \ (ri - T2) was 

neglected since (t2 — ti) ^ {ti+ T2)'^ . Noting, further, 
that if Ti -I- r2 < the pole is located out of the integra- 
tion interval, and calculating the corresponding Cauchy 
integral over the (t2 — Ti)-variable for ti + T2 > 0, one 
obtains: jfcsT'aljiV^ (^^)%2«77nF-27rru„_ A similar 

expression can be derived by expanding near the sym- 
metric point Ti 2'K + Ti and T2 — > T2 with |ri| <C 1 and 
T2 ^ 1. Therefore, the quadratic term is written as 

(Ift) u rr.^2 ^r^^^^ Ao ^ ^ 



ksTajfN- 



inp \huJc 



Ree 



Since near the poles u, ^ — 3> 1, the spatial integral 
in EqjT] is dominated by very small distances, a result 
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consistent with the locahty of the quadratic term. Also 
note that the final expression does not depend on a^. , a 
result consistent with the fact that the structure of the 
vortex lattice does not influence the quadratic term. 



B. The quartic term 

1. Useful analytical expressions 

The calculation of the next order term, the quar- 
tic term, is much more complicated since, unlike the 
quadratic term, it is strongly affected by the coupling 
of the electrons to the vortex lattice. However, the use of 
the representation, EqEl for the single electron Green's 
functions facilitates greatly the entire 8-fold spatial inte- 
gration by transforming the corresponding integrand into 
a multiple Gaussian form. Following the derivation de- 
scribed in detail in Appendix A the quartic term can be 
written as a 4D 'temporal' integral: 

r— />oo 

r24^ = nf^L^Ny^\dctM~^\^ dTidT2dT3dTi (9) 
2 Jo 

where 



1 



aiQ;2Q!3Q;4 /3o det U 



St 



(10) 

(11) 



where the vector L is given in Eg lASl and the matrix U 
in EqIMl 

The calculation of det U can be done by using the 
relations Hi + j = and noting that it can be factor- 
ized to: det [/ = Xb^a where Xa = (A3 — A2), Ab = 
^(Aa-fA^) with /3 = i(-i- + J_ + J_ + J_), A3 = 

A2 



aiQ2tt3 a2Ct3Ct4 aia20t4 aiasctji' aiQ3 ct2CK4 

Substituting these values to $ [r] one arrives at the com- 
pact expression: 



$1 



2^^ 



1 



ai 



+ q;2 + 03 + Q;4 (1 — 7^)^^^ 



(12) 



with 



_ A2 _ a2a4 — cna3 
^ A3 ai + a2 + + 



(13) 



The calculation of the exponential term results in 



1 1 

4 



-1l 


Aa ^ A 


Ab 


29o7 


{ 

V7 + I 



(14) 



Substituting Eqs. [12] and [14] to Eq. [9] wc obtain the 
final result 



h = 



1 



-ksTalN 



Ao 



hjjj, 

dTldT2dT3dT4 



(15) 



g-ro„(Ti+T2+T3+r4)-iTl_F(Ti-r2-|-T3-r4) ^ 



/3(7) 



ai -I- a2 + as -I- 04 



where: 

/?(7) = 



(1-7^)^/^ 



(16) 



^exp 

St 



O Q ^ 1 2 /^l - 7 2 1 +7,2 

—ziast On s H 1 

4^"Vl+7 1-7 



2. Major analytical properties 

The structure function /3 (7), expressed in EqfTOl con- 
trols the coupling between the four electrons involved and 
the vortex lattice. Its most remarkable feature is associ- 
ated with the dual singular points at 7 — >■ ±1, where the 
lattice sums over s or t (depending on whether 1 — 7^-0 
or 1 -I- 7 — >■ , respectively) can be replaced by inte- 
grals (over s = Sa/ (1 — 7) /2 or t = ty^ (1 -I- 7) /2, respec- 
tively) , enhancing the singularities of the corresponding 
pre-exponential factors to simple poles: 



/3(7^ 

^exp 

t 

1 



1 



1 2 



ax (1 - 7) 
1 



dse i'o'*" X 



1-7 



1-7 



1/2 



(1-7)' 
/3(7^-l) 



^exp 



1 2 



1 

ax (1+7) 
1 



~ 1 „2T2 



1 + 7 



1 + 7 



1/2 



(1 + 7) 

Note that in the sum over the remaining variable, t 
or s, only the single term t (or s) = survives, due to 
the large negative real-part values of the corresponding 
exponent. 

It is interesting to note that the values of the individ- 
ual electronic "time" variables, Tj, satisfying the singular 
conditions, 7 ±1, are given, respectively, by: 



Ti = r3 — > 0, T2 — > mr ~ t,t^ ^ nir + r 



(17) 
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T, T2 = Ti 







(18) 



where r is an arbitrary real number in the interval: 
— TT < r < TT, and n = 0,±1,±2,.... Thus, the electrons 
at such highly correlated pairs of cyclotron orbits are res- 
onantly coupled to the entire vortex lattice, yielding only 
purely harmonic contributions to the SC free energy in 
the dHvA frequency F = npH since under these condi- 
tions: e"™^(^i-^2+^-'"^*) ^-1-K%nnp _ j^Q^g g^jgQ i-j^^i- .^^ 

the singular points the factor e-'^''('^i+'^2-i-T3-i-T4) equal 
to e"^'^!"™"!, determining the thermal damping of the 
quantum oscillatory part of the SC free energy, and a 
natural (thermal) cutoff for the integrals over Tj. 

Another type of singularities of the integrand in EgfTSl 
corresponds to the vanishing denominator ai + a2 + as + 
a4 , which takes place at simultaneous vanishing of all 



a," 



1 



namely when Tj — >■ 2?™, 



0,1,2, 



At the corresponding poles the effective coupling param- 
eter of the electrons to the vortex lattice 7 — >■ , and one 
recovers the well known local approximation in which the 
electrons are only weakly coupled to the vortex lattice. 



3. Effect of the vortex lattice 

The effect of the vortex lattice on the free energy can be 
expressed more clearly by transforming the lattice double 
sum in Eall6l into a 2D reciprocal vortex lattice summa- 
tion. To do so the summation over t is transformed by 
means of Poisson formula into: 



E 



dtcxp 



2i (7rm - 



9s) t 



1-7 



1/2 




— — Trm -Os) — 

1 + 7 / V TT 



(19) 



Now, using two primitive vectors spanning the vortex 
lattice: a = xux , b = xbx + yby , with by = n/ax , the 
corresponding primitive vectors spanning the reciprocal 
vortex lattice are: a* = xby — yb^^ b* = ya-c, so that 



TT 



,f (^)' = (sa*+mb*)^ 



and 

/3(7) 



(1 + 7) 



s,m— — 00 

5a* + mb* 



1-7 
1+7 



(20) 



A similar procedure in which Poisson formula is used 
with respect to the summation over s leads to an ex- 
pression identical to Eg 1201 after exchanging 7 < — > —7. 
Since under the integrations over all , j = 1, 4 , 
a given value of 7 (corresponding to a given selection 
of Ti , T3 , r2 , r4 ) always appears with its opposite value 
—7 (corresponding to the values of ri , T3 exchanged with 
those of r2, T4) one can always replace j3 (7) in Eq llSl with 
the symmetric expression: 



G 

iv/— ^ ^^'^ 
24- 1 (1 + 



(21) 



(1+7) 



exp 



1 + 7 



(1-7) 



exp 



1 -7 



without altering the result of 14. 

Near the singular points 7 = ±1 Eql2T] describes two 
additive coherent processes of two electron pairs moving 
in cyclotron orbits on the Fermi surface and undergo- 
ing scatterings by the vortex lattice. Near the singular 
point 7 — )■ 1, where the positions of the electrons la- 
beled (1,3) along their cyclotron orbit coincide (i.e. for 
Ti = T'3 — > 0, see EaiT7|. the other two electrons, labeled 
(2,4), are moving coherently along their cyclotron orbits 
in opposite directions (i.e. T2 — > n-K — t,t^ — >■ riTr + r, 
sec Eq fT5|) . Thus, the singular 7 — > 1 contribution of 
the first term within the brackets in Eqsl21] is associ- 
ated with the electrons labeled (1,3) and involves many 
G-vectors, whereas the singular 7 — )■ 1 contribution of 
the second term is associated with the other two elec- 
trons labeled (2,4) and involves only the G = channel. 
Similarly, near the dual singular point 7 — > — 1, where 
the (1, 3) electrons arc moving in opposite directions (i.e. 
Ti rni — T. ^ niT + T ) and the positions of the (2, 4) 
electrons along their orbit coincide (i.e. T2 = T4 — )• ), 
the contribution of the first term involves only the G = 
channel, whereas the contribution of the second term in- 
volves many G-vectors. 

The physical meaning of the singular 7 — > ±1 contri- 
butions is therefore apparent: The two electrons whose 
positions on the cyclotron orbit coincide at the singular 
point undergo local mutual scattering and so exchange 
many G-vectors through the vortex lattice during the 
scattering process, while those electrons moving coher- 
ently on a large cyclotron orbit in opposite directions are 
mutually scattered through the entire vortex lattice, and 
so do not exchange momentum. 

The resulting leading contributions to /3*''™ (7) can be 
therefore written in terms of very simple formulas: The 
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forward scattering contribution takes the form: 



^0^0(7) 



1 

2 V 1 



1 



1 



1-7 



(22) 



whereas the rest of the reciprocal lattice 
contributions, ^Q_^Q /3q™ (7), which involve increasingly 
large numbers of reciprocal lattice vectors as 7 — >■ ±1, 
can be well approximated in these limiting cases by the 
two-dimensional integral: / /3q"" (7) d?G, yielding: 



of the vortex lattice on the fermionic quasi-particles at 
high magnetic field. 



4. Numerical calculations 

For numerical calculations we use Ea llSI assuming a 
square vortex lattice with ax = \/t^- Performing Poisson 
summation over s or t in Eq ll6l one can transform /3 (7) 
into simpler, equivalent forms: 



E 



G#0 



1 

1+7 



7 



1 

-1 



7 ^ 

(23) 

Note that the G = term, given by Eq. [22l represents 
the effect of the spatially uniform component of the SC 
order parameter on the free energy whereas the rest of 
the terms in Eq. [^T] correspond to all possible umklapp 
(coherent) scattering processes by the vortex lattice. 

To gain further insight into this remarkable 
coupling to the vortex lattice we may expand 
/3 (7) / («! + q;2 + as + 04) about 
gular points, say 7 — 1. and 
integrations to derive a more 
approximate) expression for 1^. 
plicity, on the first harmonic of the dHvA frequency 
F = npH, our small expansion parameters are (see 
also Sec.B4): = i (n + T2 + + r4) - 7r/2, 6 = 
\ (n - r2 + Ta - T4) + TT, ^3 = Ti - Ta, and ^4 = T4 - T2 , 
so that to second order, the key composite variables are 
given by: 



one of the sin- 
carry out the Tj- 
transparent (but 
Focusing, for sim- 



1-7 
1-^7 



1 



7*6 



(24) 



and: 



^ / d6 exp <j i(,2 



-(2C1-I-6) 

- (26-12) 



\\G\'-2np 



lr,\ 1 



Considering the umklapp scattering terms with large 
vectors G it is clear that the dominant contribu- 
tions originate from reciprocal lattice vectors satisfying: 
i |G| « \J2nF., namely having length close to the Fermi 



surface diameter. Furthermore, due to the large values of 
np and the discrete nature of G (which are measured in 
units of the magnetic length) with an elementary unit of 
about TT , the integration over ^2 yields erratically oscil- 
lating function of n^?, which reflects dramatic influence 



f^s, (7) 




1 — 7 ^ ' 



(25) 



1 + 7 
■ 7 

1-7 

1-7 
■ 7 

1 + 7 



The integrals in Eq llSl can be more conveniently eval- 
uated by transforming to the new variables (shifted 
with respect to ^i, defined above Eq [M)) : ^1 = 

^ (n + T2 + T3 + T4), ^2 = ^ (n - T2 + Ta - T4), ^3 = 
[t\ - ^-3), ^4 ^ (-T2 +T4): 
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-251 

dCs / d^i 



2?i-h«2 

X 

-(2Ci+e2) 
Hi) 



(26) 
(27) 



(25i-?2) «! + Q!2 + as + Q!4 



The distribution function I3 (^1), has been calculated 
numerically for different integer values oi np- Selecting 
integer values oi np pins the SC free energy at max- 
ima of its magnetic quantum oscillations, allowing to de- 
termine their amplitude for any given harmonic in the 
dHvA frequency F = Hnp. The result for (^1) is 
shown in Fig.l. It appears as a series of sharp peaks 
located around the points = |/c with fc = 0, 1, .., hav- 
ing monotonically increasing intensity with increasing or- 
der k. The maximum positions of the peaks are slightly 
shifted with respect to ^fc toward larger values due to the 
^i-dependence of the ^3, ^4-integrals. The peaks' height 
is found to increase with increasing harmonic order k as 
k^ , but the number of significantly contributing peaks is 
limited by the thermal damping factor e"^™''^! . A sim- 
ple estimation shows that for 2'K^kBT > AfiiOc the con- 
tribution of the second harmonic docs not exceed 10% of 
the first harmonic where the 3rd harmonic contribution 
is less than 1%. On the other hand at temperatures as 
low as 2TT^kBT < fiwc many harmonics provide compa- 
rable contributions. In this low temperature limit, re- 
placing summation over harmonics with integration one 
finds for an integer np: kBTJ2u 1^ d^ie-^"^"^^ I3 (^i) ^ 
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FIG. 1: The distribution function h (^i) given by Eal27lfor 
np ~ 32. The peak near = is too small to be observable 
in the considered scale . 



{2v+l)k 



huJr 



Therefore, the quartic term diverges as ^ as T — > 0, due 
to the resonance pairing conditions^ characterizing the 
zero spin splitting situation considered here. Note, how- 
ever, that the resulting divergence is weaker than that 
obtained in the local approximation^. In the latter the 
quartic term for an integer np was found to be propor- 
tional to ksT^^q^, with qi, 

and X, = {2u + 1 



peraturc limit: ksTj^ull ^ W' 
quadratic term, which is local in nature, is characterized 
by the low temperature limit: kBT^^q^, 



coshX„-|-cos 2ir(nF + l/2) ' 

with the following low tem- 
Note also that the 



T- 



For the sake of illustrating the novel (quantum) fea- 
tures of the theory we will focus here on the leading mag- 
netic quantum oscillatory effect by considering the first 
harmonic of the thermodynamic potential in the dHvA 
frequency F = Hnp. This situation corresponds to the 
usual dHvA experimental conditions when higher har- 
monics are relatively small. Under these circumstances 
the main contribution to I4 (see Eg 1^5)) originates in the 
second peak at ^1 ~ 7r/2, which is dominated by the in- 
tegral over small intervals around ^2 — i^r, and to lesser 
extent by all other values of ^2- The resulting integral 
over ^2 (with the integrand including e"^*"''?^ , see Ea l?7|) 
in the small intervals near ^2 — iTJ" yield the dominant 
contribution to the first harmonic. As usual for the first 
harmonic one may restrict the thermal Matsubara sum- 
mation to the single term = 0. 

The integration over ^1 around the point ^1 ~ tt/2 
has been performed for different integer values of np un- 



der the assumption that e ^^•^^^ 



-27rro„=o, The result 



presented in Fig. 2 (blue line) shows clearly the erratic os- 
cillatory dependence on np associated with the coupling 
to the vortex lattice. 

The mean base fine of this function, shown in the fig- 
ure, corresponds to I4 (^1 ~ 7r/2) calculated after replac- 



FIG. 2: The four-fold r- integral, I4,, calculated at integer 
values of np for a square vortex lattice (blue solid line). The 
red dashed line presents the harmonic part, obtained by using 
Eq. 1281 Use of the first term in this expression (corresponding 
to the single pole at 7 = 1) in the calculation of I4, yields the 
green dotted line. 



ing Psq (7) with: 



P^n^ (7) ^ 



1 - 7 1 -t- 7' 



(28) 



The result is purely harmonic, as can be seen 
by expanding one of the denominators, e.g. 
(ai 4- Q!2 + as + Q!4) (1 — 7) — >• near ti,T3 
0;T2,r4 — > TT , in the small variables (x; <^ 1), 

Tl = Tl , T2 — > — TT + T2 , r3 — > , T4 — > —1T + T4 , Up tO SCC- 

ond order and keeping only leading terms in each variable 
(see Eq l24|) . The pole contribution at ^2 (^^2 = — tt + £,2 

yields the first harmonic e2™f"^e"2"^(^i+3^') , which is 
strongly localized around the origin along both direc- 
tions ^1 and ^3 with a characteristic width ^ ySF' 
integral over ^4 is not local and it is restricted only by 
its integration limits ± (2^i — ^2) — ±27r. The remaining 
local (Gaussian) behavior in the corresponding 2D 
subspace enables one to estimate the global dependence 
of 4^'^' (a ^ 7r/2) on np as 4^"^ (^i ^ 7r/2) ^ JL. 

Fig. 2 also confirms the conclusion drawn in Sec.IIIBS 
on the basis of an analytical consideration saying that 
umklapp scattering of electron pairs by the vortex lattice 
via large reciprocal lattice vectors across the entire fermi 
sphere diameter leads to erratic oscillatory dependence 
of the thermodynamic potential on np = about the 
base line envelope ^ The absence of similar Umk- 

lapp scattering effects in the leading, quadratic term in 
the order parameter expansion, and their expected in- 
creasingly enhanced appearances in higher order terms 
of this expansion indicate that the irregularity discussed 
above should appear pronounced far from the SC tran- 
sition where the quartic and higher order terms become 
important. 

The final result for the first harmonic of the SC ther- 
modynamic potential, up to fourth order, can be written 
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FIG. 3: The "erratic" function w (np) showing the relative 
contribution to I4 associated with the coherent scattering by 
the vortex lattice. Note that the negative jumps are due 
to paramagnetic distortions of the cyclotron orbits traversing 
through vortex core regions, whereas the positive jumps are 
associated with diamagnetic distortions. 

in the form: 
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FIG. 4: 74, calculated at integer values of np for a square 
vortex lattice, as in Fig. 2, in which the singularities at 7 = 
±1 are removed (see text) with the regularization parameter 
a — .01 (red solid line). The (blue) dashed line, which is 
proportional to n.p^'''^, is the best fitting curve, whereas the 
(green) dotted curve follows a fitting formula proportional to 



r3/2 



1 - 



(29) 



1 , ^^ '^^ 
+ 0^0 (1 + W (nFj) 

2 np 



Ao 



where wq ~ 1.1 arises from the spatially uniform com- 
ponent of the SC order parameter, and is purely har- 
monic, whereas w (np), shown in Fig. 3, represents effects 
of umklapp scattering by the vortex lattice leading to de- 
viations from the purely harmonic Fourier spectrum. 

It is interesting to note that these "erratic" umklapp 
scattering processes can be viewed in real space as aris- 
ing from the passages of paired electrons in cyclotron 
orbits (near the fermi energy) through vortex core re- 
gions, where the cyclotron orbit is strongly distorted by 
the pair-potential into small circular orbits around the 
vortex cores^. The resulting deviations from the normal 
state cyclotron orbit in a vortex core are paramagnetic 
or diamagnetic, depending on the electron energy relative 
to the Fermi surface, with the paramagnetic sectors lead- 
ing to the sharp drops of the free energy shown in Fig. 3, 
while the diamagnetic ones yielding the sharp rises seen 
there. 

The existence of these erratic oscillations is due to the 
highly coherent cyclotron motions of the two pairs of 
electrons responsible for the singular terms 7 = ±1 in 
EqEH A scattering process of these electrons which can 
destroy this coherence should lead to removal of the sin- 
gular behavior. Leaving to future publications the ques- 
tion of how such scattering processes can be implemented 
into the present theory (see the discussion in Sec.V), it 
is desirable to investigate the robustness of the quartic 
term I4 with respect to smearing of the singularities at 
7 = ±1. This can be done by artificially shifting 1 — 7, 
and 1 + 7 in Ea l25l slightly away form their vanishing 
forms to 1 + (T — 7, and 1 + cr + 7 respectively, for small 



values of (T > 0, and repeating the calculation shown in 
Fig. 2. The result for a = 0.01 is shown in Fig. 4. In 
addition to the significant reduction of the overall mag- 
nitude and suppression of the (coherent-scattering) "er- 
ratic" oscillations, the np dependence of the mean base 



line changes from np to n 



-3/2 



characterizing the local 



approximation of the GGL theoryi^. 



IV. THE EFFECT OF VORTEX LATTICE 
DISORDER IN THE WHITE NOISE LIMIT 



Calculation of the influence of vortex-lattice disorder 
on the SC free energy in the magneto-quantum oscilla- 
tions limit can be performed analytically in the white 
noise limit. Invoking the general expansion of the state 
function ipo{x, y) in terms of Landau orbitals wave func- 
tions, ip[){x, y) = e"^ Yl,n c„e*'^"^~^^+*"/^) , the structure 
factor takes the form: 



Hi) 



1 



1 



(1 _ ^2)1/2 
exp 



(30) 



E 

nst 



1+7 



1+7^2 



where the coefficients {c„} may be considered as random 
variables. Averaging the structure factor over realiza- 
tions of these coefficients and exploiting the usual (Wick) 
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decoupling: 

(/3(7)> = 



a. (1 - 72)1/2 



E 



exp 



7 2 

-s 



1 + 7 



1 + 7 ^2 

1-7 ■ 



, in the white noise hmit, i.e. (c*c„+s) (c*^_s_,_tC„+t) - 
<^ri+.s+t,n+s'5n,n+t = <5t,o , One finds: 



(/3(7)) = 



(1-7^) 



1/2 



(31) 

The final step in the procedure leading to the white 
noise limit should be the replacement of the discrete 
summations in Eg 1311 with integrations (e.g. by taking 
Oa; — > oo there), resulting in the expression: 



(/3(7)) 



1 



1 



(32) 



1—7 1+7 

This is a rather surprising results since it is seen to 
be twice Eql22| obtained for the forward scattering term. 
The latter (i.e. the G = term), which is usually associ- 
ated with all incoherent scattering processes, is expected 
to be the sole survivor of an averaging over white-noise 
disorder, and as such to coincide with Eql32l 

In this limiting case, only incoherent scattering pro- 
cesses by the vortex matter contribute to the SC ther- 
modynamic potential, and the final result, up to fourth 
order, is purely harmonic, with the first harmonic given 
by: 



i]ii''))Mi'^)^i-. 



r3/2 



An 



1 TT-" 

2 np 



Ao 



(33) 

i.e., very close to the well known Maki-Stephen 
expressioni^ii^, as expanded to the same order in Aq. 

An interesting question arises here as to wether the 
white-noise average of higher order terms in the order- 
parameter expansion presented in this paper also agree 
with the self-consistent Born approximation (SCBA) in- 
herent to the Maki-Stephen approach^. In particular, 
possible destruction of the highly coherent motions of 
the electron pairs responsible for the singular contribu- 
tions to the quartic term by an infinite subset of di- 
agrams which are topologically distinct from the quartic 
diagram, might lead to significant deviations from the 
SCBA. 



Landau theory of strongly type superconductivity at high 
magnetic fields. An exact analytical expression for the 
quartic term in the corresponding order parameter ex- 
pansion, having a physically transparent form, is pre- 
sented. The resulting expression reveals singular non- 
local contributions to the SC thermodynamic potential, 
associated with highly coherent cyclotron motions of 
the paired electrons near the Fermi surface, which are 
strongly coupled to the vortex lattice. The dominant 
contributions to the SC free energy, arise from incoher- 
ent scattering by the spatially averaged pair-potential, 
which is purely harmonic in the dHvA frequency. How- 
ever, coherent scatterings by the ordered vortex lattice 
generate, at low temperatures, erratically oscillating (i.e. 
paramagnetic-diamagnetic) contribution to the SC free 
energy as a function of the magnetic field, associated with 
sharp distortions of the large quasi-particle cyclotron or- 
bits on the Fermi surface traversing through vortex core 
regions. Vortex lattice disorder, which tends to suppress 
this oscillatory component, is found to simplify consid- 
erably the calculation allowing analytical evaluation of 
higher order terms in the order-parameter expansion. 
However, it can be shown that the infinite subset of di- 
agrams constituting the standard, self consistent Born 
approximation (SCBA)^, exploited in the white noise 
limit of the discorded vortex system^, have the same 
type of singular points as that found in our calculation 
of the quartic term. It would be therefore very interest- 
ing to search for, and then evaluate subsets of diagrams, 
topologically distinct from those appearing in the SCBA, 
which might, after resummation, destroy the highly co- 
herent cyclotron motions responsible for the above singu- 
larities. Physically speaking, it is expected that the effect 
of impurity-scattering on the paired electrons, as calcu- 
lated beyond the relaxation time approximation, could 
destroy this coherence. Whether or not the robustness 
of this type of singularities with respect to scattering of 
quasi particles by a disordered vortex matter is destroyed 
by going beyond the SCBA is a crucial question in our 
understanding of the vortex lattice disorder on the dHvA 
oscillations in the SC statei^. 
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CONCLUSION AND DISCUSSION 



Appendix A 



A novel Green's function representation is exploited in 
this paper for a microscopic derivation of the Ginzburg- 



Similar to the calculation of the quadratic term, it is 
convenient to introduce the following, center of mass and 
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relative coordinates: 



R = 




+ r2 


+ r3 


+ r4) 


Q = 




r2 


+ r3 


-1-4) 


D = 




r2 


r3 


+ r4) 


P = 




+ r2 


rs 


-1-4) 



where pi = Ti — ri_i. This transformation can be 
written in the matrix form: X = M * r, where X = 
{R, Q, D, P} are four 2D vectors and A/ is a 4 x 4 ma- 
trix with |det A'/| = 1/2. 



= 7 (P1-P2+P3-P4) 



J 



An ingredients of the quartic term, which depend on the electronic spatial coordinates, i.e.: 

^4 = rif / d2{r}f4({r})i^4({r}). 



n 



(0) 



2tt 1 



{2n) 



Ao 



Ki{{Y}) = / dTidT2dr3dT4e-*(^i-^^+^='-^*)"^-^-(^i+^^+^^+^'') 
Jo 



1 



■ CXp [- {p.ipl + P2pI + /i3P3 + A*4pI)] 



aia2a3a4 

r4({r}) = .g*(ri,r2).g(r2,r3)g*(r3,r4).g(r4,ri)</Jo(ri)<^S(r2)</'o(i-3)</'S(r4) 

will be rewritten now in terms of the new coordinates. Let us start with the gauge factors, 
S*(ri,r2)3(r2,r3)g*(r3,r4)5(r4,ri) = e''", where: 



'7s 2 '^['"1 ^ "^2] - [r2 X rg] + [rg x r4] - [r4 x ri]) 2i {Q^Ry - QyRx) , 

which depends only on the vectors R, Q. 

The product of the four Landau orbitals, labeled by ni = n + s + t, ?i2 = ?^ + s, 713 = n , and 714 = n + i , is given 
by the following expression: 



(m) (sg) (!m) (0) 

?7a +'?a +'?a 



with: 



(m) 
»?A 



??A 



„(0) 
'^A 



i ^ Ej-Xjl/j = i [2Q^Ry + 2QyR^ + {D^Py + P^Dy)] 

3 

-T.y' = -{Ql + ^Rl + {Dl + p';)) 

3 

3 

* + ('3-^ ^4 - '^D^t - 2P^s) - go {RyNi - DyS - P^^t) 



1 2 \ ' 2 

-4102^^3 

3 



ql [2n^ + 2ns + 2nt + st + + t 



2 I .21 



and: N4 = J2 "-j = 4n + 2.s + 2t. 

The last factor, Ki{{r} ^ exp (?7g), arising from the transitional invariant parts of the Green functions, is indepen- 
dent of the center of mass coordinates: 



(mi +f^2+ /i3 + Pi) Q^ + '2 {P2 - Pi) {DQ) + 
2 {pi - p3) (PQ) + {p2 + Pi) + {pi + p-i) P2 



(Al) 
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The simplest integration to carry out, over , yields a non-vanishing result only if ^ {ejiij) = 0, justifying the 
parametrization of Uj chosen above. In this case J dR^ = = cl^N^. Next, the i?y-integration. 



J exp ['iiQa:Ry - 4i?^ - qoNRy] dRy ^^y/^cxp 



(A2) 



leads to a space independent correction, j^QqN'^, which removes the n-dependence of ry^'': 77^' + j^QqN'^ = 
— jQq [s^ + 1^)- As a result, summation over Landau orbitals is trivially done, yielding the total number of orbitals 

^n^ = Ny. 

Combining the Q-dcpendcnt terms the corresponding integral is: 

/ d^Q exp [-/3o [Q^ + 2/324 {DQ) + 2/3,3 (PQ)]] - ^ cxp [/3o {Pl^D^ + Pf^P^ + 2^24^13 (DP))] (A3) 
J Po 



where 



« 1 , , 1 1 « M2 - ^ (mi - Ms) 

Po = Ml + M2 + M3 + M4 + 1, P24 = ^ , Pl3 - 



/?0 



/So 



(A4) 



The most complicated analytical part of the calculation, the £)P-integrations, is now done by introducing the 4D 
vectors: 



Z = {D^,Dy,Px,Py} ,L = qo{-it,s,-is,t} 



and the r dependent 4x4 matrix: 



(A5) 



U = 



(M2 + M4) - 



-/3/324/3l3 

-i/2 



-/3/?24/3i3 -V2 

(M2 + M4) - /3/3I4 + 1 -V2 -/3/?24/3l3 

-^/2 + - 

-/3/324/3i3 (mi + Ms) - /3/3?3 + 1 



(A6) 



and then performing the resulting Gaussian integrations to have: 



/ d^Zexp [-Z'^UZ + LZ] = exp [i^[/-^L] 

J V dct U 



(A7) 
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